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Abstract 

Some aspects of applications of bunching parameters are dis- 
cussed. It is investigated to what extent Monte-Carlo models, which 
have been tuned to reproduce global event-shape variables and 
single-particle inclusive distributions, agree with each other. 
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1 Introduction 

One of the simplest observable, which contains information about the dy- 
namics of multiparticle production beyond single-particle densities, is the 
multiplicity distribution. While the study of the multiplicity distribution 
P/v m full phase space deals with limited dynamical information influenced 
by charge- and energy-momentum conservation, the investigation of the 
evolution of the probabilities P n (5) of detecting n particles in ever smaller 
sizes S of phase-space windows (bins) can provide detailed information on 
QCD multihadron production without these trivial constraints. A devia- 
tion of this distribution from that expected for purely independent particle 
production can be attributed to dynamical local multiplicity fluctuations. 

The important quest behind such a study is the understanding of 
the origin of short-range correlations between final-state particles, leading 
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to the appearance of dynamical multiparticle spikes in individual events. 
As a consequence of these correlations, the normalized factorial moments 
(NFMs) 

oo 

{n [q] ) = Y,n [q] P n (5), n [q] = n(n - 1) . . . (n - g + 1), (2) 

n=q 

of the local multiplicity distribution P n (5) exhibit a power-like increase 
with decreasing 5, namely F q (5) oc 5^^ q Q. The constants g are called 
intermittency indices. This phenomenon reflects the peculiarity of P n (5) 
to become broader with decreasing 5. Since NFMs satisfy the scaling 
property F q (XS) = \^^ q F q (5), this is widely regarded as evidence that the 
correlations exhibit a self-similar underlying dynamics. 

Experimentally, local fluctuations in e + e~-processes have already been 
studied by the TASSO, HRS, CELLO, OPAL, ALEPH, DELPHI and L3 
Collaborations || . The data do exhibit approximate power-like rise of the 
NFMs with a saturation at small 5. The conclusion has been reached that 
such a phenomenon is a consequence of the multi-jet structure of events, 
i.e., groups of particles with similar angles resulting in spikes of particles 
as seen in selected phase-space projections. The hard gluon radiation sig- 
nificantly affects the NFMs, so that they have stronger increase in 3-jet 
events than in 2-jet events. It has been found that for the statistics used at 
that time current Monte-Carlo models can, in general, describe the data, 
even without additional tuning. 

Recently, it has been realized that the factorial-moment method poorly 
reflects the information content of local fluctuations, since the NFM of 
order q contains a trivial contamination from lower-order correlation func- 
tions (see reviews ||). As a result, rather different event samples can 
exhibit a very similar behavior of the NFMs. The fact that subtle details 
in the behavior of P n (5) are missing, together with the small statistics 
used, may be the reason why different Monte-Carlo models can reasonably 
describe the local fluctuations measured in e + e~ annihilation so far. 

Another shortcoming of the factorial-moment measurement is that in 
moving to ever smaller phase-space bins, the statistical bias due to a finite 
event sample (N ev ^ oo) becomes significant, especially for high-order 
moments q. This is because in actual measurements the NFMs at small 
bin size are determined by the first few terms in (0). In most cases this 
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leads to a significant underestimate of the measured NFMs with respect 
to their true values. 

Cumulants are a more sensitive statistical tool (see and references 
therein). However, their measurement is rather difficult and was rarely 
attempted. Besides, the cumulants are expected to be influenced by the 
statistical bias to even larger degree, since they are constructed from the 
factorial moments of different orders q. 



2 Local Properties 

An important step towards an improvement of experimental measurements 
of the local multiplicity distribution was made in 0,0, where it was shown 
that any complex distribution can be represented as 

\n n 

P n {5) = P (5)-L n , L n = Uvr +1 (5), 

where A = Pi(5)/Po(5). The factor L n measures a deviation of the distri- 
bution from a Poisson with L n = 1. Non-poissonian fluctuations exhibit 
themself as a deviation of L n from unity. The L n is constructed from the 
bunching parameters (BPs) 

q P q (5)P q - 2 (5) 
^ = ~q^l PUV ' q {) 

The values of the BPs and NFMs for most popular distributions are shown 
in Table 1. The most interesting observation is that while the NFM is an 
"integral" characteristic of the P n (S) and the BP is a "differential", both 
tools have values larger than unity if the distribution is broader than a 
Poisson. Generally, however, one should not expect that all BPs are larger 
than unity for a broad distribution; BPs probe the distribution locally, 
i.e. they are simply determined by the second-order derivative from the 
logarithm of P n (5) with respect to the n. Note, that in the case of local 
distributions, the width of distributions is mainly determined by r]2(S). 
This observation is based on the simple fact that P n {5) ceases to be bell- 
shaped at sufficiently small 5. 
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Distribution 


Pn 


NFMs 


BPs 


Pos. Binomial 

Poisson 

Neg. Binomial 

Geometric 


C*p n (l - P ) N ~ n 
p n exp(—p)/n\ 

r(n+k) n(i | \-(k+n) 

r(n+i)r(fe)^ v ^ p > 
p n (p+l)- n ~ 1 


n? =1 (i-i)<i 
i 

m = i(i+i)>i 
n? =1 (i+*) > i 


q-2-N ^ 
1 

q-l+k ^ -I 
q-2+k ^ 

9-1 



Table 1. NFMs and BPs for positive-binomial, Poisson, negative-binomial 
and geometric distributions. 

BPs are more sensitive to the variation in the shape of P n (S) with 
decreasing 5 than are the NFMs ||. In the case of intermittent fluctua- 
tions, one should expect rj 2 (5) oc 5~ d2 . For multifractal local fluctuations, 
the r) q (5) are ^-dependent functions for all q > 3, while for monofractal 
behavior r] q (5) = const for q > 3 [Q. 

From an experimental point of view, the BPs have the following im- 
portant advantages [Q: 

1) They are less severely affected by the bias from finite statistics than 
the NFMs, since the gth-order BP resolves only the behavior of the multi- 
plicity distribution near multiplicity n = q — 1; 

2) For the calculation of the BP of order q, one needs to know only the 
g-particle resolution of the detector, not any higher-order resolution. 

The problem we are dealing with in this paper is to investigate whether 
different Monte-Carlo (MC) models, which were tuned to reproduce the 
global-shape variables and single-particle inclusive densities, can lead to 
the same structure of the local multiplicity fluctuations which are deter- 
mined by many-particle inclusive densities. We study JETSET 7.4 PS 0, 
ARIADNE 4.08 § and HERWIG 5.9 [|] models. The models have been 
tuned as described above by the L3 Collaboration fTCf] . 
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3 Monte- Carlo Analysis 



1) Horizontal BPs: 

In order to reduce the statistical error on the observed local quantities 
when analyzing experimental data, we use the bin-averaged BPs @,[5|]: 

9K 1 q-l iV ? 2 _i(M) ' V ' 

1 M 

N q (M) = -Y,N q (m,6), (5) 

m=l 

where N q (m, S) is the number of events having q particles in bin m and M = 
A/5 is the total number of bins (A represents the size of full phase-space 
volume). To be able to study non-flat distributions, like for rapidity, we 
have to carry out a transformation from the original phase-space variable to 
one in which the underlying density is approximately uniform, as suggested 
by Bialas, Gadzinski and Ochs ||11||. 



2) Generalized integral BPs: 

To study the distribution for spikes, we will consider the generalized 
integral BPs || using the squared pairwise four-momentum difference 
Q\ 2 = — (pi — V2) 2 ■ I n this variable, the definition of the BPs is given 
by 

Xq{Q) ~q-l II^CQS) ' W 

where U q (Q 2 ) represents the number of events having q spikes of size 
Q 2 in the phase-space of variable Q\ 2 1 irrespective of how many par- 
ticles are inside each spike. To define the spike size, we shall use the 
so-called Grassberger-Hentschel-Procaccia counting topology for which a 
many-particle hyper-tube is assigned a size Q 2 corresponding to the maxi- 
mum of all pairwise distances (see for details). For purely independent 
particle production, with the multiplicity distribution characterized by a 
Poissonian law, the BPs @ are equal to unity for all q. 

3.1 In rapidity variable 

In order to study fluctuations inside jets, in most investigations the fluc- 
tuations have been measured in the rapidity y defined with respect to the 
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thrust or sphericity axis 0]. The Monte Carlo analysis for this variable is 
performed in the full rapidity range | Y |< 5. Fig. 1 shows the results for 
the BPs (||) for rapidity after the Bialas-Gazdzicki-Ochs transformation. 
The second-order BP for JETSET model decreases with increasing M up 
to M ~ 20, which is found to correspond to the value of M at which the 
maximum of the multiplicity distribution P n {8) first occurs at n — 0. At 
large M, all BPs show a power- law increase with increasing M, r] q ~ M aq . 
This indicates that the fluctuations in y defined with respect to the thrust 
axis are multifractal scale invariant. 
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Figure 1. BPs as a function of the number M of bins in rapidity defined 
with respect to the thrust axis. The shaded areas represent the statistical 
and systematical errors on the JETSET predictions. 

Note that the conclusion that fluctuations have a multifractal structure 
is possible without the necessity of calculating the intermittency indices 
<p q . In contrast, to reveal multifractality with the help of the NFMs, one 
first needs to carry out fits of the NFMs by a power law. 





6 



HERWIG predictions (dashed lines) significantly overestimate the second- 
order BP obtained from LUND MCs. Since, for small phase-space cells, the 
second-order BP is determined by the dispersion of the distribution ||,|5]], 
this means that the HERWIG produces too broad local multiplicity distri- 
butions. Such a result confirms that obtained by the ALEPH Collabora- 
tion [0. 

To study the disagreement between Monte-Carlo models in more detail, 
one can split 7/2 into two BPs: 

(±±) . (+-) 

V2 = m +V2 ■ (7) 

Here 772 is defined by (|J) with A^(?n, 8) = (m,8y), N2 (m,8y) 

being the number of events having like-charged two-particle combinations 
inside bin m of size 8y. Analogously, r]2 + is constructed from the number 
of events (m, 8y) having unlike-charged two-particle combinations. 

Note that due to a combinatorial reason, < rf^ \ 




Figure 2. The second- order BP as a function of the number M of bins in 
rapidity defined with respect to the thrust axis for like-charged and unlike- 
charged particle combinations. 

Fig. 2 shows that r)^^ and t]2 + indeed behave completely differently. 
While ^2 ±± ' ) shows the expected rise, r^ -1-- ' 1 shows a strong decrease at low 
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M and the onset of an increase at large M. The structure of 772 observed 
in Fig. 2 is a combination of these two effects. 

Note that 772 is strongly effected by the Bose-Einstein (BE) interference 
incorporated into the JETSET generator^. This is not unexpected since 
r]2 ~ P2/P1, which is very similar to the correlation functions used for 
Bose-Einstein studies. 

Let us remind that in order to model the BE interference in JETSET, 
the momenta of identical final-state particles are shifted to reproduce the 
expected two-particle correlation function. The main disadvantage of such 
ad hoc method is that it spoils overall energy-momentum conservation 
and it is necessary to modify also momenta of non-identical particles to 
compensate for this. This effect in JETSET model can be seen in Fig. 2. 

The strong anti-bunching tendency seen for unlike- charged particles at 
M < 30 can be attributed to resonance decays and to chain-like parti- 
cle production along the thrust expected from the QCD-string 
model fl3| . The latter effect leads to local charge conservation with an al- 
ternating charge structure. Evidence for this effect was recently observed 



by DELPHI [14|. As a result, there is a smaller rapidity separation be- 



tween unlike-charged particles than between like-charged and r]^~ ^ is much 
larger than 772^ a ^ small M. Having correlation lengths 5y ~ 0.5 — 1.0 in 
rapidity, the resonance and the charge-ordering effects, however, become 
smaller with increasing M. 

Note that to distinguish the NFMs calculated for different charge com- 
binations in a bin-splitting technique is difficult due to insufficient sensi- 
tivity of this tool and a purely combinatorial reason. 



3.2 In the four- momentum difference 

The study of BPs described above can help us to understand a tendency of 
the particles to be grouped into spikes inside small phase-space intervals. 
Another question is how the multiplicity of these spikes fluctuates from 
event to event when the spike size goes to zero. To study this, we will use 
the BPs defined in ([|). 

Fig. 3 shows the behavior of Xq as a function of — InQ 2 . The full lines 
represent the behavior of the BPs in the Poissonian case. In contrast, all 
BPs obtained from the Monte Carlo models rise with increasing — In Q 2 

1 Here and below, we show JETSET predictions with the BE interference disabled 
after the retuning of this model to describe global-shape variables. 
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Figure 4. Generalized second-order BP as a function of the squared four- 
momentum difference Q 2 between two charged particles. 
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(decreasing Q 2 ). This corresponds to a strong bunching effect of all or- 
ders, as expected for multifractal fluctuations. The anti-bunching effect 
(x q < 1) for small — \nQ 2 is caused by the energy-momentum conserva- 
tion constraint ||. 

To learn more about the mechanism of multiparticle fluctuations in 
Q12 variable, we present in Fig. 4 the behavior of the second-order BP as 
a function of — In Q 2 for multiparticle hyper-tubes (spikes) made of like- 
charged and those of unlike-charged particles, separately. A significant 
difference is observed for like-charged combinations between HERWIG and 
LUND MCs. 



4 Discussion 

Local multiplicity fluctuations in Monte Carlo models have been studied by 
means of bunching parameters. Since all high-order BPs show a power-like 
rise with decreasing the size of phase-space interval, none of the conven- 
tional multiplicity distributions given in Table 1 can describe the local 
fluctuations observed in the MC models. 

For e + e~ interactions, one can be confident that, at least on the parton 
level of this reaction, perturbative QCD can give a hint for the under- 
standing of the problem. Analytical calculations based on the DLLA of 
perturbative QCD show that the multiplicity distribution of partons in 
ever smaller opening angles is inherently multifractal Qualitatively, 
this is consistent with our results on the BPs for rapidity. Quantitatively, 
however, the QCD predictions disagree with the e + e~ data and MC mod- 



els [ 16 



In this paper we show that the power-law behavior of BPs is mainly due 
to like-charged particles. JETSET gives the same power-law trend even 
without the BE effect. This means that the intermittency observed for 
like-charged particles appears to be largely a consequence of QCD parton 
showers and hadronization. 

The predictions of the ARIADNE 4.08 model are comparable with those 
of the JETSET 7.4 PS model. This is essentially due to the same im- 
plementation of hadronization, which is based for both models on string 
fragmentation. 

A noticeable disagreement, however, is found between LUND and HER- 
WIG models. The conversion of the partons into hadrons in the first mod- 
els is based on the Lund String Model [13]. However, the hadronization 
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in HERWIG is modelled with a cluster mechanism ||. This can be a 
rather natural candidate to explain the observed difference between local 
fluctuations in these models. A particular concern is the large difference 
between MC's for r] 2 . The behavior of r)2 for not very small intervals is 
sensitive to low-multiplicity events, for which hadronization details could 
play a significant role. 
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